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Abstract 

A concept of neighborhood in complex networks is addressed based on the cri- 
terion of the minimal number os steps to reach other vertices. This amounts to, 
starting from a given network Ri, generating a family of networks Ri,l = 2,3, ... 
such that, the vertices that are i steps apart in the original are only 1 step apart 
in Rn. The higher order networks are generated using Boolean operations among the 
adjacency matrices Mi that represent Ri. The families originated by the well known 
linear and the Erdos-Renyi networks are found to be invariant, in the sense that 
the spectra of Mi are the same, up to finite size effects. A further family originated 
from small world network is identified. 



Several properties of complex networks have been addressed recently, much 
inspired by the identification of their relevance in the description of the re- 
lations among the individual constituents of many systems, which have their 
origin from natural to social sciences, telephone and internet to energy sup- 
ply, etc. [1,2,3,4]. The huge development in this research area came along with 
the proposition of many parameters that might be relevant to the character- 
ization of properties of networks. The sofar most established and quantifiers 
are the distribution of links n{k)., clustering coefficient C, mean minimal dis- 
tance among the nodes < d >, diameter D, the assortativity degree [5,6]. The 
evaluation of these and other indices for actual networks help characterizing 
them, putting into some network classes with well defined properties, such as 
the small- world [7], the scale- free [8], the Erdos-Renyi [9] or random networks, 
etc. 

As the number of nodes directly connected to node i is ki, n{k) characterizes 
the immediate neighborhood of the network nodes. In this work we explore 
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further neighborhood properties of networks, which are related to the distri- 
bution of the number of second, third, .... neighbors. For the sake of simplicity, 
we assume that the networks we consider herein are connected, i.e., each node 
can be reached from any other one. Two nodes are i — th neighbors, or neigh- 
bors of order 0{£), when the minimal path connecting them has i steps. Then, 
for a given network R, the explicit evaluation of the distributions of neighbors 
0{i) along the network, £ = 1,2, 3, ... promptly indicates the structure of min- 
imal paths connecting the nodes. This classifies uniquely the neighborhood of 
a vertex, in the sense that if two vertices are neighbors, they are not 
0{j) neighbors for £ ^ j. Also, we consider that any vertex is 0(0) neighbor 
of itself. 

It is expected that the neighborhood properties change with £. However, if a 
meaningful criterion G can be devised, according to which the ^-neighborhoods 
of it! remain invariant, it may be important to assign R a, G neighborhood in- 
variant (NI) property. It characterizes an invariance of R with respect to length 
scale, provided length is measured by number of steps along the network. Also, 
it is distinct both from the scale free distribution of vertex connections, as well 
as from geometrical scale invariance included in the construction of the net- 
work as, for instance, in the class of Apollonian networks (AN). 

In our investigation we identified, for each £, all 0{£) neighbors of R and 
constructed a family of networks R. Each i?£ G R is defined by the same set 
of nodes, while links are inserted between all pairs of vertices that are 0{£) 
neighbors in R = Ri. Thus, the family R characterizes the neighborhood 
structure of R. The Rf 's were actually set up by a corresponding family M of 
adjacency matrices (AM) Mi [5], achieved by the systematic use of Boolean 
{B) operations [10] among matrices. 

To define a neighborhood invariance criterion, it is necessary to identify a 
relevant network property that may be present in all elements of R, e.g., 
the corresponding values for ni{k) and Ce [11]. However, these are global 
indices that do not provide a sufficiently precise identification of a network. A 
much more precise characterization is based on the eigenvalue spectra of the 
family M^[12,13], even if it is well known that there may exist topologically 
distinct networks which share the same spectrum (co-spectral), and that only 
the complete set of eigenvalues and eigenvectors univocally characterizes a 
network. A spectrum based invariance criterion condenses (in the sense of 
reducing from N"^ elements of M to iV numbers) a lot of information about 
the network, and has been used herein. This can be further justified by the 
fact that, if two AM's in have the same spectrum, they differ only by a 
basis transformation. 

Before proceeding with details of the guidelines of our actual work, we recall 
that the first matrix Mi, which describes the original network has only O's 
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or I's as entries for all of its elements. If Mi is applied to a unitary vector ve, 
with all but one entry set to 0, the resulting vector expresses which vertices 
are linked to the the vertex i. If we take the usual matrix product of Mi 
by itself, the non-zero elements (Mf)^^ of the resulting matrix M^ indicates 
how many possible two-step walks along the network, with endpoints r and s 
exist. Contrary to what happens with Mi. Mf has many elements {Mf)rs > 1, 
indicating multiplicity of paths starting at r and ending at s. In particular, all 
elements {Mf)rr of the diagonal can have this property, since they count all 
two-step walks that start at r, visit any of the vertices s to which r is hnked 
{Mi)rs — 1, and turn back to r. The same interpretation is valid for all usual 
powers of Mi. 

As the elements of Mi are all O's or I's, we can regard (Mi),.^ as B variables, 
and use the B sum, subtraction and product operations [10], respectively 

e,e,®, 

0©0=01©0-10©1=11©1=1 

o©o = oi©o = ioei = oiei = o (i) 

0®0=0 100=0001=0 



to define B operations between matrices of B elements. The B matrix opera- 
tions are defined by using the usual matrix element operation rules, replacing 
the usual sum, subtraction and product involving matrix elements by the 
corresponding B operations. To avoid multiplicity of notation, we will use 
hereafter the same symbols ©, © and to indicate matrix B operations. 

If we consider M2 = Mi Mi and compare it to Mf , we see that the posi- 
tion of all their zero elements coincides, while if we collapse to 1 all non-zero 
elements of Mf we obtain M2. In fact, M2 indicates the possibility of two- 
step walks, while it deletes the information on the multiplicity of walks. As 
the neighborhood concept does not take path multiplicity into account, the B 
operations are well suited to define the matrices Mj. 

For instance, M2 can be expressed by 

M2 = (Ml © M2) - (/ © Ml) = (7 © Ml) Ml - (7 © Ml), (2) 



where 7 indicates the identity matrix. To see this note that all forward- 
backward walks, included together with pairs of distinct sites linked by two- 
step walks in M2, can not be present in M2, as any vertex has been defined to 
be 0(0) neighbor of itself. Thus we must subtract 7 from M2. Also, it is neces- 
sary to B sum Ml to M2, and subsequently subtract Mi, as M2 may describe 
two-step walks between two sites that were already related in the original net- 
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work Ri. Noting that Mq = I, Eq.(2) can be generalized for arbitrary value 
of e by: 

£—1 (.—1 l—l 

Mi = (0 Mj) <8, Ml - (0 Mj) = M,_i <g, Ml e (0 Mj). (3) 

j=0 j=0 j=0 

Once a precise procedures to set up all M^'s is available, let us briefly comment 
on several possibilities opened by the knowledge of M for the purpose of 
evaluation of network indices. 1) If we have a finite network with N vertices, 
then there is a large enough ij^ax such that M^ = 0, > imax- Thus, the value 
for D is found when the first M^ = is found. Also, when i approaches D, the 
M^'s become sparser and, as a consequence, the number of zero eigenvalues 
increases largely. 2) As for each r, s pair, {mi)rs = 1 for only one M^ e M, 
one can collapse in a single matrix 

^max 

M=Y: jM, (4) 

j=0 



all information on the neighborhood of any pairs of vertices. Particularly, all 
pairs r, s of the 0{i) neighbors satisfy M^s = 3) To obtain the average 
minimal path for each node r, it is sufficient to sum all elements of the r — th 
row (or column) of M and divide by N — 1 . The average minimal path for the 
network follows immediately. 4) The evaluation of all by means of line to 
line multiplication of matrix elements, can be easily computed once the family 
M has been obtained. 5) The matrix M can be easily used to visualize the 
structure of network with the help of a color or grey code plots. 

The numerical evaluation of the eigenvalue spectra A^, r G [1, N] for the family 
Me has been carried out for several networks. For each one of them, we are 
particularly interested to understand how the form of the spectral density 
Pe{\) as function of A evolves with i. We first illustrate the procedure by 
showing how the spectra of the standard AN depends on i. Several properties 
of AN's have been discussed recently [14,15]: they are constructed according to 
a precise geometrical iteration rule that grant them topological self-similarity 
and a scale free distribution of nodes degrees. The Mi spectra of the AN for g 
successive generations converge very quickly to a well defined form. However, 
for the finite size networks up to (? = 8, with 3283 vertices, we can not identify 
any sort of invariance in the eigenvalues density distribution ^^(A). This is 
better visualized with the help of the integrated spectra [16] 

1 } 

n.(A) = - J MX')dx', (5) 

— oo 
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Fig. 1. Plots of n X A for Mi (solid) and M2 (dashed) of the AN, when g = 7 and 8 

for successive generations, as shown in the Figure 1. Hi and 112 converge very 
quickly to distinct g independent forms. Thus, scale free distribution of node 
degrees and self similarity do not necessarily lead to NI. 

Exact NI invariance is observed for the linear chain network with periodic 
boundary conditions, where each vertex interacts only with its two next neigh- 
bors. The corresponding AM Mi has a well known pattern, where all 1 ele- 
ments are placed along the nearest upper and lower diagonals to the matrix 
main diagonal. This matrix has been used to describe a very large number 
of models in 1 dimension, like the system of phonons, tight-binding electrons 
[17,18], etc. Also, it has been used as the starting point to construct small 
world networks, by changing some of the original nearest neighbors links ac- 
cording to a given rewiring probability Pw For Mi, p is expressed analytically 
by the relation [17,18,19] 

P^^^ ~ (4-P)V2 - (6) 

The other matrices of the family Mi keep essentially the same shape, with 
two sequences of I's along near-diagonals. Each one of them moves one step 
away from the main diagonal upon increasing the value of £ by 1, but this 
does not change the spectrum. Indeed, this operation can be regarded, e.g. 
in the analysis of tight-binding systems, as a decimation procedure of half 
of the sites along with a renormalization of the hopping integral [19]. The 
resulting system has exactly the same shape as the original and, consequently, 
the same spectral density. The diameter of the network depends linearly on 
N, as diagonals move away one step when i is increased by 1. 

The Erdos-Renyi networks [9] constitute another class where one could expect 
to find NI. Indeed, if connections are randomly distributed for Ri, so should 
they also be for all members of R. For very low values of p (the connection 
probability between any two nodes), Ri is split into several disjoint clusters, 
and this situation is preserved for all other Rgs. In such cases, the Mi spectral 
density does not obey a simple analytical expression, being constituted by 
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Fig. 2. Spectral density pi for the Erdos-Renyi network single realizations, (a) 
N = 1000,p = 0.05, £ = 1,2,3,4,5. (b) Scaled p2 for conditions p = cjN^ when 
c = 68, z = 0.9, = 500, 1000, 3000, 3500 
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Fig. 3. Plots of the spectra p x A of the matrices Mi^l = 1,3,5,7,9, for a 
Watts-Strogatz small world network starting from a nearest neighbor linear chain 
and pyj = 0.02. 

some individual peaks superimposed on a shallow background. The M( (i >2) 
spectra also share the same qualitative structure. Quantitatively, it is observed 
an increase of the A = dominant peak and a decrease of the other peaks. 
This is related to the clustered structure of the network and to the fact that 
several clusters are reaching their own diameter. 

When p = c/N^, with z < 1, and for such large enough c, almost all nodes 
are connected in a single cluster, and p obeys the well known semicircle law 
[20,21] 

p(A) ~ (4 - XY'. (7) 



In this regime, networks usually have a very small diameter, and invariance 
can only be noted for a few values of i. For the average node number (k) ~ 0.5, 
we have found that the p2 also obeys (7), as shown in the Fig. 2b. However, 
for smaller values of (k), a clear skewness in the distribution is observed for 
P2, despite the semicircle form of pi. We conclude that Erdos-Renyi networks 
are NI for a restricted subset of p and c. 
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Fig. 4. Almost invariant spectra p x A of the matriees M£,l = 7—10, for a 
Watts-Strogatz- small network starting from a next-nearest neighbor linear chain 
and py, = 0.2. 

We have also investigated small world networks obtained along the Watts- 
Strogatz rewiring procedure [7], starting both from the above discussed next- 
neighbor (nn) and the nearest-next-neighbor (nnn) linear chains. In Figure 3, 
we show a sequence of Mi spectra for very small value of p^, starting from a 
nn chain. We see that the p(A) is split into two parts: the first one, essentially 
described by Eq.(6), corresponds to the contribution of unperturbed segments 
of the linear chain. We see that, as it happens to the whole spectra when 
Pw — 0, this part remains almost invariant, at least for a large number of M^. 
The other part of the spectra, located outside of this first region, has no well 
defined shape. It entails a considerably lower number of eigenvalues, which 
increases with i. This indicates that, at each increase of i, a small number 
of eigenvalues migrates from the unperturbed part into it, as the number of 
vertices that are affected by the rewiring operation increases with i. This 
successive migration ends up by affecting the whole form of the spectrum for 
large enough i. Of course this behavior depends on the value of p^. If we 
increase it, there is a smooth transition on the shape of the Mi spectrum, 
until it reaches a pattern with many structures and bands, similar to that of 
the AN. 

For p^ in the range p^ ~ 0.15 — 0.22, we have observed that, as £ increases, the 
spectral density pi evolves towards a very peculiar form, which remains almost 
invariant for several values of i, as shown in the Figure 4. This spectrum has 
its own features, distinct from those discussed before for the fully ordered 
and disordered networks. For the specific = 0.21 and N = 1000, we 

have a large diameter D = 15, and the shown form of the spectra remains 
invariant for i e [7, 10]. For smaller values of i, the shape changes steadily 
from a structured shapes similar to those in Figure 3 into the invariant form. 
For larger values of i, finite size effects lead to quite sparse M^, with a large 
number of zero eigenvalues: pe evolves to a 5 like distribution centered at A = 0. 

This effect can be graphically illustrated with the help of the matrix M. In 
Figure 5 we draw the position of the 0{£) neighbors for three distinct ranges of 
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Fig. 5. For the same Watts-Strogatz small network of Figure 4, graphical illustration 
of distribution of O(^) neighborhoods: I € [1,6], [7, 10], [11, 15] for (a), (b) and (c) 
respectively. 



i. The particular shape in Figure 4 is associated with roughly dense matrices 
when i e [7,10]. For smaller and larger values of matrices have rather 
distinct structure. 

For other values of N, we have observed the same evolution. For instance, 
when N = 1500 and 2000, the shape lasts almost invariant for more gener- 
ations, respectively i e [8, 12] and [11, 16], indicating that this behavior can 
be more robust as N increases. For p^^, larger than the range given above, this 
persistence in the form is not observed. For smaller values of p,^, the spectra 
changes very slowly as shown in Figure 3. In such cases, finite size effects set 
in prior than any tendency of evolution towards the form shown in Figure 4. 

To conclude, in this work we have discussed the concept of higher order neigh- 
borhood and neighborhood invariance of networks. These have been obtained 
by a systematic use of Boolean matrix operations and the definition of a AM 
family. We explored well-know networks, showing that this property is not 
equivalent to other concepts of scale and geometrical invariance. Further, we 
looked for evidence of NI based on the invariance of the spectral density, iden- 
tifying this property in the linear chain, Erdos-Renyi network, and finally, in 
a non-trivial class that evolves from Watts small world network. 
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